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The ground state of F has been observed as an unbound resonance 220 keV above the ground state of ^^F. 
Comparison of this result with USDA/USDB shell model predictions leads to the conclusion that the F ground 
state is primarily dominated by .sii-shell configurations. Here we present a detailed report on the experiment in 
which the ground state resonance of " F was first observed. Additionally, we report the first observation of a 
neutron-unbound excited state in F at an excitation energy of 2500(220) keV. 
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I. INTRODUCTION 

The neutron-rich region around N = 20 has been a topic of 
active experimental and theoretical research for over 30 years, 
owing to the transition from pure sd to mixed sd-pf shell con- 
figurations first deduced from mass measurements of neutron- 
rich sodium isotopes [1, 2]. As the available intensities at rare- 
isotope beam facilities have increased, it has become possible 
to explore increasingly neutron-rich systems near A^ = 20, in- 
cluding those at and beyond the neutron dripline. 

The heavy fluorine isotopes represent some of the most 
neutron-rich N ^ 20 systems that can be measured with 
present experimental techniques. A 2004 measurement [3] 
reported two 7-ray transitions in each of the 25,26,27p jgotopes, 
with the higher-lying ^^F transition confirmed in a recent ex- 
periment [4]. In ^^F, a transition was observed at 777(19) keV 
and assigned to the first 1/2+ excited state. This is in poor 
agreement with USD [5, 6] shell model predictions which 
place the 1/2+ at 1997 keV. SDPF-M Monte Cario Shell 
Model calculations [7], which allow for sd-pf shell mixing, 
are in better agreement with observation, placing the 1 /2+ at 
1 100 keV. This suggests that the first 1/2+ excited state in ^^F 
exhibits significant sd-pf configuration mixing. Additionally, 
Ref. [3] reports a low-energy transition not predicted by USD 
in each of the 25.26,27p isotopes, speculating that these tran- 
sitions might correspond to 1/2^ states arising from proton 
p-sd cross shell excitations. 

Until recently only one measurement of neutron-unbound 
states in fluorine isotopes has been reported. A 28 keV reso- 
nant decay from ^^F was assigned to a 1/2^ excited state in 
^^F at an excitation energy of 4249 keV [8]. In the present 
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paper we report the observation of the unbound ground state 
of ^^F and an unbound excited state in ^^F. The results of the 
^^F experiment have been reported in a recent article [9]. 



II. EXPERIMENTAL DETAILS 
A. Setup 

The experiment was performed at the National 
Superconducting Cyclotron Laboratory (NSCL) at Michigan 
State University. Unbound states in 27,28p y^gj-g populated 
by nucleon removal from a beam of ^^Ne. The ^^Ne beam 
was produced by first accelerating '*^Ca^''+ to 140 MeV/u in 
the NSCL coupled cyclotrons [10]. The "^^Ca then impinged 
upon a 1316 mg/cm^ ^Be production target. Products of the 
calcium on beryllium reaction were sent through the A 1900 
fragment separator [11], which was tuned to optimize the 
transmission of ^^Ne at 62 MeV/u. The A 1900 included an 
achromatic aluminum wedge at its second image point to 
disperse fragments according to A/Z and improve separation. 

After the A 1900, a quadrupole triplet magnet focused the 
beam onto a 288 mg/cm^ ^Be reaction target. Upstream of 
the target, the beam passed through a pair of position-sensitive 
cathode readout drift chambers (CRDCs) separated by 227 cm 
and a pair of plastic scintillators separated by 1044 cm. The 
location of each of these detectors along the beam axis is 
shown in Fig. 1 . The CRDC position measurements were used 
to calculate the beam position on the reaction target by ray 
tracing through the quadrupole triplet. The upstream scintil- 
lator was 1010 /xm thick, and the downstream ("target") scin- 
tillator was 254 jim. Each scintillator recorded a time sig- 
nal, and these signals were used to calculate the beam time of 
flight (th). Additionally, the target scintillator recorded an en- 
ergy loss signal (AE^,). As shown in Fig. 2, the vaiious beam 
components were well separated in energy loss versus time of 
flight. The desired ^^Ne composed approximately 2% of the 
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FIG. 1. (color online) Diagram of the experimental setup. 
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FIG. 2. Energy loss versus flight time of the secondary beam, show- 
ing its three components: ^^Ne, ■'^Mg, and various lighter species 
("wedge fragments") produced in the A1900 wedge. 



put was indicative of the fragment energy. 

Due to the size and complexity of the setup, separate 
data acquisition (DAQ) systems were used for MoNA and 
Sweeper-CAESAR. Events in each DAQ were recorded with 
a timestamp, allowing coincidences to be reconstructed off- 
line. Although run separately, the triggering of each DAQ was 
controlled by a shared logic module, which allowed for trig- 
ger conditions involving both subsystems. To reduce dead- 
time, the experiment required coincidences between MoNA 
and the 5 mm scintillator located at the back of the Sweeper 
box. CAESAR detected y rays in coincidence, but they were 
not a required trigger condition. 



B. Data Analysis 



1. Charged Particle Separation 



beam, and the remainder was composed of ^^Mg (87%) and 
various lighter species produced in the aluminum wedge. 

A diagram of the experimental setup is shown in Fig. 1. 
The experiment consisted of three subsystems, each used to 
measure a different type of reaction residue potentially re- 
sulting from the breakup of neutron-unbound states in 27:28p. 
neutrons, y rays (from feeding to bound excited states in the 
daughter), and residual charged paiticles. Neutrons were de- 
tected in the Modular Neutron AiTay (MoNA) [12], which 
measured their time of flight, position, and the amount of light 
deposited, y Rays were detected in the Caesium Iodide Array 
(CAESAR) [13], which measured their total energy and time 
of flight. Charged particles were first deflected 43° by the 
Sweeper magnet [14]. They were then detected in a pair of 
CRDCs, an ionization chamber measuring energy loss, and 
two plastic scintillators. The front face of each scintillator 
was 40 cm x 40 cm, and each was coupled to four photo-tubes. 
The upstream scintillator was 5 mm thick and recorded a time 
signal; this signal was combined with the time output of the 
target scintillator to determine the fragment time of flight. The 
downstream scintillator was 150 mm thick, and its charge out- 



The charged particle measurements allowed for event-by- 
event isotope identification after making a variety of correc- 
tions to the data. The first step was to identify the various el- 
ements reaching the end of the Sweeper using measurements 
of energy loss and total energy. Energy loss (AEf) was ob- 
tained from the ionization chamber signal. The fragment time 
of flight {tf) and charge output of the 150 mm scintillator (Q) 
were each used as an independent indicator of total energy. 
Figures 3(a) and 3(b) show the element separation in AEf-Q 
and AEf-tf, respectively. In the final analysis, events were 
required to fulfill conditions in both parameter spaces. 

For a given element, isotopes were separated by construct- 
ing a corrected time of flight parameter (f(.) indicative of A/Z. 
The corrections to the time of flight accounted for the vary- 
ing paths taken through the Sweeper, and the primary indica- 
tors of this path length were the dispersive position (x) and 
dispersive angle {Ox) of the fragment as it exited the magnet. 
Additionally, a variety of other parameters (c.f. Table I) were 
found to correlate with the time of flight of a given isotope 
and were included in the corrections. Due to the lack of focus- 
ing elements, as well as non-homogeneities in the Sweeper's 
magnetic field, it was necessary to consider three-dimensional 



TABLE I. Correction factors used for isotope separation. To calcu- 
late the corrected time of flight, we take the sum of each factor multi- 
plied by its corresponding parameter and then add this sum to tf. The 
symbols x{y) and 6x{dy) respectively refer to the dispersive (non- 
dispersive) position and angle of the charged fragment as it exits the 
Sweeper. The symbol xtrgtiytrgt) denotes the beam's dispersive (non- 
dispersive) position on the reaction target. The remaining symbols 
are all introduced in the text. 



Parameter 



Correction Factor 



X (mm) 

x^ (mm^) 

X (mm ) 

9r (mrad) 

e'} (mrad^) 

9;? (mrad^) 

xOx (mm mrad) 

xQx (mm mrad^) 

X 0x (mm mrad) 

x^Qx (mm^ mrad^) 

y^ (mm^) 

6y (mrad) 

xtrgt (mm) 

ytrgt (mm) 

tb (ns) 

Q (arb.) 

A£ (arb.) 



-5.0595 X 10- 

-8.97 X lO""* 

-3.0x10"^ 

+8.0x10^2 

-1.0x10^5 

+2.0x10 

-1.5x10 

-6.0x10 

-2.0x10 

+ 1.4x10 

+ 1.0x10 

-3.0x10 

+1.7x10^2 

+4.0 X 10"^ 

+1.0x10^' 

+1.3 xlO"^ 



-6 

-4 
6 
6 

7 
3 
3 



h4.0 X 10 



-3 



correlations between tj, x, and 0,, as well as non-linearities, in 
determining the appropriate corrections. Because of superior 
statistics, the time of flight corrections were determined for 
fluorine elements produced from the ^^Mg beam. These same 
corrections were then used to separate the isotopes of interest, 
fluorines produced from ^^Ne. 

As shown in Fig. 3(c), a three-dimensional plot of tf-x-Ox 
displays isotope bands. For the purpose of time of flight cor- 
rections, it is useful to reduce the x-9x phase space into a sin- 
gle "emittance" parameter e(x, Ox), as this will allow for cor- 
rections to the flight time to be made in a straightforward way. 
To determine e{x,9x), the tf-x-Qx scatter-graph was profiled 
by dividing the x-Ox phase space into small regular rectangu- 
lar regions and finding the mean tf for each region [15]. This 
profile plot is shown in Fig. 3(d), with the grayscale level rep- 
resenting mean tf. From here, the location of Ox as a function 
of X was fit along the lines of constant tf in the profile. As 
shown by the curve in Fig. 3(d), the location of these lines 
was well-described by a second order polynomial. 



f{x)= ax + bx + c [mrad] , 



(1) 



with a = 0.010391 mrad/mm- and b = 0.84215 mrad/mm. 
The final constant c can take on any value; it only causes the 
curve to shift to a different line of constant tf. 

Once f{x) was determined, e{x, Ox) was constructed simply 



As shown in Fig. 3(e), plotting e{x, Ox) versus tf reveals iso- 
tope bands in two dimensions. From here, an initial corrected 
time of flight parameter was calculated by projecting onto the 
axis perpendicular to the bands. The time of flight corrections 
were then further refined by iteratively removing any correl a- 
tions between tf and the parameters Usted in Table I. 

The final corrected time of flight (f^ ) for fluorines produced 
from ^^Ne is shown in Fig. 3(f). By fitting this spectrum with 
the sum of five Gaussians constrained to have equal width, 
we determined the 26p.27p cross-contamination to be approxi- 
mately 4%. The factors used in constructing the corrected time 
of flight are listed in Table I, and it should be noted that the 
most important corrections (in addition to those for x. Ox, and 
their higher order combinations) are those for y^ and x,rgt. 



2. Decay Energy Calculation 

The decay energy of the breakup of unbound states was 
calculated using invariant mass analysis. In Euclidian coor- 
dinates, the decay energy E^ is expressed as 



E,= 



- 2 [EfEn - PfPn COS 0) - nif - 



(3) 



eix,Ox)^Ox-fix) 



(2) 



where mflm,,), Ef{E„), and Pf{pn) refer to the mass, energy, 
and momentum of the charged fragment (neutron), respec- 
tively, and is the opening angle between the two decay prod- 
ucts. Charged fragment inputs to Eq. 3 were determined us- 
ing a partially inverted COSY transformation matrix [16, 17], 
which operated on the measured position and angle behind the 
Sweeper and the x position of the beam on target. The trans- 
formation returned the energy and angle at the reaction target, 
as well as the track length and the target y position. 

The neutron input to Eq. 3 was calculated from time of 
flight and position measurements in MoNA using relativis- 
tic kinematics. The trigger logic was designed such that the 
stop for each MoNA time digitizer channel was provided by a 
delayed signal from the target scintillator Thus the recorded 
time signals were a measurement of neutron time of flight (r„). 
To calibrate the raw digitizer signals, a linear slope and offset 
were applied to each channel. The slopes were determined 
from a pulser run; relative offsets between MoNA bars were 
determined from cosmic-ray muon tracks; and an overall off- 
set was set from the travel time of prompt y rays. Vertical 
and lateral positions in MoNA were assumed to be at the cen- 
ter of the interaction bar, and the horizontal position was cal- 
culated from the time difference between signals measured 
on either end of the bar In the case of multiple interactions 
within MoNA, the earliest hit with f„ > 40 ns was used in the 
analysis. The cutoff of 40 ns was chosen to eliminate any ran- 
dom first hits that arrived too early to be prompt neutrons. 

A plot of the neutron time of flight to the front face of 
MoNA is presented in Fig. 4, for three conditions: ungated 
(including all incoming beam components), ^^F produced 
from ^^Ne, and ^^F produced from ^^Ne. In the ungated plot, 
the peaks from prompt neutrons and y rays are clearly identi- 
fiable on top of a random flat background consisting primar- 
ily of room background yrays and cosmic-ray muons. When 
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FIG. 3. a) Element separation in AE-E. b) Element separation in AE-ti,. Panels (a) and (b) are both gated on incoming Ne and have the 
fluorine events circled, c) Three-dimensional plot of /^ vs. dispersive position and angle after the Sweeper for fluorine fragments produced from 
the ■'^Mg beam. The various bands correspond to different isotopes of fluorine, d) Profile of the three-dimensional plot in (c), including the 
emittance parameter, e{x, 0,), that is determined by fitting lines of constant f/, with a second order polynomial, e) Plot of the e(x, Ox) parameter 
determined in (c) versus f/,, demonstrating isotope bands in two dimensions, f) Final corrected time of flight {tc) for fluorine elements produced 
from the Ne beam, with F and F indicated. 



requiring coincidences with 26-27p^ jj^g jj^^ background is es- 
sentially eliminated, and the prompt neutron peak dominates 
the spectrum. 



3. y-Ray Measurements 

CAESAR was calibrated using a variety of standard 7-ray 
sources [18]. Although a large magnetic shield was placed 
between it and the Sweeper, CAESAR was still subject to sig- 
nificant fringe fields (on the order of 3 mT) which affected the 
response of its photo-tubes. To account for this, the array was 
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FIG. 4. Neutron time of flight to tiie front face of MoNA, for tiie 
following conditions: a) ungated; b) in coincidence with F; c) in 
coincidence with F. 



calibrated with the Sweeper set to the rigidity at which the ex- 
periment was performed. Furthermore, to account for poten- 
tial hysteresis effects, a recalibration run using a '^^Y source 
was taken any time the field of the Sweeper was changed dur- 
ing the experiment. 

For 7-ray events depositing light in multiple crystals, the 
total deposited energy was calculated using an "add-back" 
technique [13]. The in-beam data were Doppler corrected; 



for the correction, the detection point taken as the center of 
the first interaction crystal, and the emission point assumed 
to be the center of the reaction target. To reduce background 
from random coincidences, only events falling within a spe- 
cific time window were included in the final analysis. Because 
of electronic effects (walk in the leading-edge discriminators 
used for timing), the time window was implemented as a two- 
dimensional cut on time of flight versus Doppler-corrected en- 
ergy. 



C. Modeling and Simulation 

Resonant states were modeled by a Breit-Wigner line-shape 
with an energy dependent width derived from 7?-Matrix theory 
[19]. The equation for the line-shape is 



(7(£;£o,ro,£) 



ATt{E-To) 



[E(, + ^l{E■To)-EY + \[Tt{E■To)] 



(4) 

where A is an amplitude, Eq is the central resonance energy, 
Fq parameterizes the central resonance width, (. is the orbital 
angular momentum of the resonance, and F^ and Af are given 
by 



Yt{E)^2Pi{E)yl 
A({E)^~[S({E)-SdEQ)]yl 



(5) 



The P( and S( functions in Eq. 5 are related to the spherical 
Bessel Functions, J('{p), and their derivatives; 






(6) 



^1/2 



with Fe = {np/2y''Jt+i,2{p) and Gf = (-1)' (7rp/2) 



^1/2 



J 



(Mil) 



(P)- 



In addition to resonant states, a non-resonant background is 
expected in the ^^F -^^^ F + « decay energy spectrum, result- 
ing from the decay (via emission of a neutron with E^i < 3 
MeV) of high-lying continuum states in ^**F to high-lying 
states in ^^F that subsequently feed the ground state of ^^F. 
The ^^¥ fragment can then be detected in coincidence with 
the first neutron, giving rise to the background distribution. 
This background was modeled as a Maxwellian distribution 
of beam velocity neutrons. 



/(e;0)=AJe/03e 



-£/© 



(7) 



with the temperature a free parameter. This model provides 
a good fit to the observed non-resonant data and has been em- 
ployed in number of other invariant mass measurements, for 
example [20-24]. 

Broadening due to experimental resolution and acceptance 
was accounted for in a Monte Carlo simulation of the exper- 
iment. In the simulation, the kinetic energy of the incoming 
^^Ne beam was modeled as a Gaussian with £0 = 62. 1 MeV/u 



and Or = 1 .72 MeV/u, clipped at £ < 64.5 MeV/u. The beam 
angle and position were also modeled as Gaussian with Ox = 
11 mm, Oqx = 4.0 mrad, C7y = 9.0 mm, and Ggy =1.1 mrad. 
Additionally, the dispersive angle and position were given a 
correlation of O^/x = 0.0741 mrad/mm. The angle and po- 
sition of the incoming beam were determined from position 
measurements in the two CRDC detectors upstream of the re- 
action target. The beam energy was determined by comparing 
measured and simulated distributions in the two downstream 
CRDC detectors for runs where the reaction target was re- 
moved. The ^Be(^^Ne,^^'^^F) reactions were treated in the 
Goldhaber Model [25] including a small friction term [26] to 
degrade the beam energy by 0.6%. The transport of charged 
fragments through the Sweeper was simulated using a third 
order COSY transformation matrix, produced from measure- 
ments of the Sweeper's magnetic field [27]. 

The resolution of charged particle position and angle mea- 
surements was modeled as Gaussian, with Gpos = 1.3 mm and 
<7a«g = 0.8 mrad. These resolutions were determined from 
data taken with a tungsten mask shadowing the CRDC de- 
tectors. The primary acceptance cut concerning the charged 
particles was the requirement that they pass through the 5 = 
±150 mm active area of the downstream CRDC. Neutron time 
of flight resolution was modeled as Gaussian with a — 0.3 ns, 
and the neutron x-position resolution was modeled as a sum 
of two Laplacian functions; 



Pi-- 



2C71 



(I-Pl 



2C72 



(8) 



with ai = 16.2 cm, 02 = 2.33 cm, and pi = 53.4%. The 
form of Eq. 8 and the parameters CJi, CJ2, and pi were de- 
termined from shadow bar measurements and GEANT3 sim- 
ulations [28]. As mentioned, the neutron y and z positions 
were assumed to be at the center of the detection bar, result- 
ing in a uniform uncertainty of ±5 cm. The overall resolution 
and acceptance for the decay of ^^F into ^^F + n has already 
been presented in Ref. [9], and the corresponding shapes are 
essentially identical in the case of ^^F* breakup. 

Due to the low statistics of the present data set, an unbinned 
maximum likelihood technique was used for parameter esti- 
mation [29]. This technique involves forming a small range, 
Ri, around each experimental data point and then summing 
the number of weighted Monte Carlo points that lie within 
the volume. To marginalize systematic errors resulting from 
Monte Carlo fluctuations and non-linearities within the Ri, the 
generated model sets were made large (^3x10^ events), and 
the volume size was chosen to be small (0.05 MeV). 



III. RESULTS AND DISCUSSION 
A. 27f Excited State 

The black squares in Fig. 5 show the measured decay en- 
ergy spectrum of ^^F + n coincidences. As mentioned previ- 
ously, we expect a non-resonant contribution in the ^''F + n 
data, so they were fit with the sum of a Maxwellian distri- 
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FIG. 5. (color online) Measured relative energy spectrum for ^^F + « 
coincidences. The filled squares with error bars are the experimen- 
tal data, the dashed red curve is the result of a 380 keV resonant 
simulation, the shaded grey curve is a simulation of the Maxwellian 
non-resonant background (0 = 1.48 MeV), and the solid black curve 
is the sum of the resonant and non-resonant models, with a reso- 
nant/total fraction of 33%. The inset is a plot of the negative log- 
likelihood as a function of central decay energy, with each point min- 
imized with respect to all other free parameters. 



bution and an £ = 2 Breit-Wigner resonance, using the tech- 
nique outlined in Section II C. In the fit, the resonance energy 
£0, resonance width Fq, Maxwellian temperature 0, and res- 
onant/total fraction /, were all allowed to vary freely. In order 
to extract £0 , the parameter of interest, a profile log-likelihood 
was constructed by scanning a range of £0 values and plotting 
the negative log-likelihood (— ln[L]) minimized with respect 
to the other free parameters (Fq, 0, and /). This profile like- 
lihood curve is displayed in the inset of Fig. 5, and it reaches 
a clear minimum at £0 = 380 keV. The na confidence inter- 
vals were determined from the In [£„,„,/£] > «^/2 limits. As 
indicated on the figure, the lc7 and 2cJ confidence intervals 
were ±60 keV and jjq keV, respectively. The best-fit values 
of the other parameters were determined to be Fq = 10 keV, 
= 1.48 MeV, and / = 33%. The simulated best fit curves 
are superimposed on the data in Fig. 5, with the dashed red 
curve representing the 380 keV resonance, the shaded grey 
curve the Maxwellian background, and the solid black curve 
their sum. A comparison between simulation and data is also 
shown for selected intermediate parameters (neutron time of 
flight, fragment kinetic energy, neutron-fragment opening an- 
gle, and neutron-fragment relative velocity) in Fig. 6. 

The presumption of £ = 2 decay is based on a pure single- 
particle model in which the least-bound neutron resides in the 
0^3/2 shell. In reality, configuration mixing and shell evolu- 
tion could lead to significant contributions from decay with 
other orbital angular momenta. Separate analyses using i^ I 
and i = 3 resonances yield results that do not differ signifi- 
cantly from the i = 2 case. The lack of sensitivity to i val- 
ues is largely due to experimental resolution, which is lim- 
ited primarily by uncertainty of the reaction position within 
the ^Be target. The width of the measured resonance is al- 
most completely determined by experimental response, over- 
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FIG. 6. (color online) Comparison of ^^F + n simulation and data 
for a) neutron time of flight; b) charged fragment kinetic energy; c) 
neutron-fragment opening angle; d) neutron-fragment relative veloc- 
ity. The filled black squares are the experimental data, and the solid 
black, dashed red, and shaded grey curves depict the same simulation 
components as in Fig. 5. 
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FIG. 7. Summary of experimentally known levels in 26.27p^ includ- 
ing the present observation of an unbound excited state in F at 2500 
keV, decaying to the ground state of ^^F. The shaded grey boxes 
around the various levels indicate the total uncertainty in their place- 
ment relative to the F ground state. The dashed lines surrounding 
the presently observed 2500 keV level represent the uncertainty on 
the decay energy only, and the total uncertainty also includes that of 
the ^^F In separation energy. All bound excited information is from 
[3], ground state energies are from [30], and ground state J^ are from 
[31]. 



shadowing any differences that might arise from varying the 
I value. Contribution from £ = decays might also be possi- 
ble, but such decays cannot be separated from the Maxwellian 
background since the resolved lineshape of the two models 
is very similar for small absolute scattering lengths (|flj| < 5 
fm). A scattering state near threshold (larger \a^\) is clearly 
not present since the data display no enhancement at low de- 
cay energy. 

Only two counts were observed in CAESAR in coincidence 
with 26F + n {Ey = 760 and 1180 keV). In the case of 100% 
branching to a bound excited state in ^^F, roughly 50 counts 
would be expected in CAESAR, based on the approximate 7- 
ray detection efficiency of 30% [13]. Thus the observation 
of only two 7 rays in CAESAR indicates that the presently 
observed decays feed the ground state of ^^F, allowing for 
an unambiguous assignment of the observed resonance to an 
excited state in '^F The most recent mass measurements of 
26'27f [30] place the ^^F ground state 2120(210) keV above 
the ground state of ^^F, so we assign the presently observed 
380(60) keV resonance to a 2500(220) keV excited level in 
2'F. Fig. 7 presents this newly observed level along with the 
other measured states in 26'27p p^ 4^ -50, 31]. 

To interpret our observations, we have performed shell 
model calculations using the USDA and USDB interactions 
[32], which operate in the traditional sd model space (Oiis^, 
1.51/2, ^nd 0(i3/2 for both protons and neutrons). The calcula- 
tion results are compared with experiment in Fig. 8. As seen 
in the figure, each calculation predicts three or more states in 
the same energy region as our observation. Extending the cal- 
culations to include pf shell components would only compli- 
cate the situation since opening up the model space increases 
the available number of excited state configurations. The as- 
signment of the observed resonance to a specific state is not 
possible because the reaction {\p-\n removal) does not pref- 
erentially populate one state over the others. 
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FIG. 8. Measured excited levels in ^^F compared to predictions of 
the USDA and USDB shell models. The first two measured excited 
states are both from [3], while the state at 2.5 MeV is from the present 
work. 
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FIG. 9. (color online) Relative energy spectrum for ^^F + n coinci- 
dences. The filled squares with error bars are the experimental data, 
the dashed red curve is the the 220 keV resonance simulation, and the 
dotted blue curve is the 810 keV simulation. The solid black curve 
is the sum of the 220 keV and 810 keV resonances, with the relative 
contribution of the 220 keV resonance at 28%. The inset shows the 
profile log-likelihood as function of the lower resonance energy. 



B. 28f Binding Energy 

As discussed in Ref. [9], the measured ^^F decay energy 
is best described as a sum of two independent i = 2 Breit- 
Wigner resonances, with the lower resonance at 220(50) keV 
(Fo = 10 keV), the upper resonance at 810 keV (Fq = 100 
keV), and the lower resonance composing 28% of the total 
area. As with ^^F, the width of each measured resonance was 
dominated by experimental resolution, making sensitivity to 
the resonance £ value minimal. An ^ = scattering state was 
excluded based on incompatibility with the measured data, 
and a non-resonant Maxwellian background was not expected 
since ^^F was populated directly by one-proton knockout from 
^^Ne. The measured ^^F decay energy spectrum is presented 
in Fig. 9, along with the best fit two-resonance simulation 
and the profile log-likelihood curve. Additionally, Fig. 10 
shows a data-simulation comparison for neutron time of flight. 
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FIG. 10. (color online) Comparison of ^^F + n simulation and data 
for a) neutron time of flight; b) charged fragment kinetic energy; c) 
neutron-fragment opening angle; d) neutron-fragment relative veloc- 
ity. The filled black squares are the experimental data, and the solid 
black, dashed red, and dotted blue curves depict the same simulation 
components as in Fig. 9. 
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FIG. 1 1 . (color online) Difference between experimental and theoret- 
ical (USDA, USDB) binding energies for A' = 19 isotones, 9 < Z < 
17. The error bars on the data points represent experimental errors 
only. The blue dotted and red dashed bands represent the respective 
170 and 130 keV RMS deviations of USDA and USDB interactions. 
Experimental values, save for Z = 9 which is from the present work, 
are taken from [30] if reported there; otherwise they are from the 
2003 Atomic Mass Evaluation [33]. Figure reproduced from Ref. 
[9]. 



fragment kinetic energy, neutron-fragment opening angle, and 
neutron-fragment relative velocity. No 7 rays were recorded 
in CAESAR in coincidence with ^^F+«, and around 30 would 
be expected in the case of 100% branching to excited ^^F. This 
indicates that the observed resonances feed the ^^F ground 
state. 



The present observation of the ^^F ground state as a 
220(50) keV unbound resonance can be combined with the 
^^F mass measurement of Ref. [30] (^^F atomic mass excess 
equal to 24630(190) keV) to calculate the ^^F binding energy 
as 186040(200) keV. By comparing measured binding ener- 
gies with the predictions of the UDSA/USDB shell model, 
which does not allow for mixing between sd and pf shell con- 
figurations, it is possible to qualitatively determine the contri- 
bution of pf shell "intruder" components in the ground state 
of a given nucleus. Such a comparison is shown in Fig. 1 1 for 
A^ = 19 isotones with 9 < Z < 17. As seen in the figure, the 
agreement is very good for the heavier isotones closer to sta- 
bility (Z > 13), while it becomes dramatically worse for the 
isotones with 10 < Z < 12 which lie within the island of in- 
version. At Z = 9, the good agreement between USDA/USDB 
and experiment is dramatically recovered, indicating that in- 
truder components play a minimal role in the ground state 
structure of ^^F. This suggests the existence of a low-Z bound- 
ary (or "shore") of the island of inversion beginning at Z = 9. 



IV. CONCLUSIONS AND OUTLOOK 

In conclusion, we have used the technique of invariant mass 
spectroscopy to make the first determination of the ^^F bind- 
ing energy at 186040(200) keV. Additionally, we have ob- 
served a neutron-unbound excited state in neighboring ^^F, 
with 2500(220) keV excitation energy. 

Interpretation of the ^^F state in terms of shell model pre- 
dictions is difficult due to the large number of levels predicted 
near 2500 keV and uncertainty in the reaction mechanism 
used to populate ^^F*. The level structure of ^^F is relevant 
to a variety of open questions in nuclear physics, including 
the transition from pure sd to mixed sd-pf neutron config- 
urations and its associated consequences (such as the large 
oxygen-fluorine dripline shift of six or more neutrons [34]). 
Additionally, it has been suggested [3, 35] that proton p-sd 
cross-shell excitations could play a role in the structure of 
low-lying ^^F excited states, possibly in tandem with sd-fp 
shell breaking on the neutron side. As such, it would be inter- 
esting to revisit unbound excited states in ^^F experimentally, 
using a direct reaction mechanism that can selectively pop- 
ulate specific states. Possible reactions include one- or two- 
proton knockout (from ^^Ne or ^^Na) and ^^F(ii, p) in inverse 
kinematics. 

The measured ^^F binding energy indicates a low-Z bound- 
ary of the island of inversion at A^ = 19. It would be interesting 
to further explore this mass region to see if this trend con- 
tinues. Extension of the present technique to lighter A^ = 19 
isotones (Z < 8) would be very difficult, if not impossible, 
since they are all unbound by three or more neutrons [36- 
40] . However, a similar technique could potentially be used 
in the N = 20 isotonic chain by performing a direct mass mea- 
surement of bound ^''F. For this purpose, the precision obtain- 
able with time-of-flight techniques at current in-flight radioac- 
tive beam facilities would likely be sufficient. Such a mea- 
surement would be particularly interesting since the SDPF-M 
Monte Carlo Shell Model predicts ^^F to have a very large in- 
truder occupation of 91.5% (62.7% two-particle, two-hole ex- 
citation and 28.8% four-particle, four-hole) [41]. Measuring 
its mass would provide the first experimental data on ^^F for 
comparison with theory and help to better explain the evo- 
lution of shell structure in the low-Z (< 10) region around 
A? = 20. 
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